Friction vs Texture at the Approach of a Granular Avalanche

We perform a novel analysis of the granular texture of a granular bed close to stability limit. Our analysis is based on a unique criterion of friction mobilisation in a simulated two-dimensional packing. In this way, we recover the bimodal character of granular texture, and the coexistence of weak and strong phases in the sense of distinct contacts populations. Moreover, we show the existence of a well-defined subset of contacts within the weak contact network. These contacts are characterized by their important friction, and form a highly coherent population in terms of fabric. They play an antagonistic role with respect to force chains. We are thus able to discriminate between incoherent contacts and coherent contacts in the weak phase, and to specify the role that the latter plays in the destabilisation process.Comment: 4 pages, 6 figure

Continuum simulation of the discharge of the granular silo: a validation test for the mu(I)-visco-plastic flow law

Using both a continuum Navier-Stokes solver, with the mu(I)-flow-law implemented to model the viscous behavior, and the discrete Contact Dynamics algorithm, the discharge of granular silos is simulated in two dimensions from the early stages of the discharge until complete release of the material. In both cases, the Beverloo scaling is recovered. We first do not attempt quantitative comparison, but focus on the qualitative behavior of velocity and pressure at different locations in the flow. A good agreement is obtained in the regions of rapid flows, while areas of slow creep are not entirely captured by the continuum model. The pressure field shows a general good agreement. The evolution of the free surface implies differences, however, the bulk deformation is essentially identical in both approaches. The influence of the parameters of the mu(I)-flow-law is systematically investigated, showing the importance of the dependence on the inertial number I to achieve quantitative agreement between continuum and discrete discharge. The general ability of the continuum model to reproduce qualitatively the granular behavior is found to be very encouraging.Comment: 12 pages, 15 figure

The granular silo as a continuum plastic flow: the hour-glass vs the clepsydra

The granular silo is one of the many interesting illustrations of the thixotropic property of granular matter: a rapid flow develops at the outlet, propagating upwards through a dense shear flow while material at the bottom corners of the container remains static. For large enough outlets, the discharge flow is continuous; however, by contrast with the clepsydra for which the flow velocity depends on the height of fluid left in the container, the discharge rate of granular silos is constant. Implementing a plastic rheology in a 2D Navier-Stokes solver (following the mu(I)-rheology or a constant friction), we simulate the continuum counterpart of the granular silo. Doing so, we obtain a constant flow rate during the discharge and recover the Beverloo scaling independently of the initial filling height of the silo. We show that lowering the value of the coefficient of friction leads to a transition toward a different behavior, similar to that of a viscous fluid, and where the filling height becomes active in the discharge process. The pressure field shows that large enough values of the coefficient of friction ($\simeq$ 0.3) allow for a low-pressure cavity to form above the outlet, and can thus explain the Beverloo scaling. In conclusion, the difference between the discharge of a hourglass and a clepsydra seems to reside in the existence or not of a plastic yield stress.Comment: 6 pages, 6 figure

Multiscale Analysis of the Stress State in a Granular Slope in Transition to Failure

By means of contact dynamics simulations, we analyze the stress state in a granular bed slowly tilted towards its angle of repose. An increasingly large number of grains are overloaded in the sense that they are found to carry a stress ratio above the Coulomb yield threshold of the whole packing. Using this property, we introduce a coarse-graining length scale at which all stress ratios are below the packing yield threshold. We show that this length increases with the slope angle and jumps to a length comparable to the depth of the granular bed at an angle below the angle of repose. This transition coincides with the onset of dilatation in the packing. We map this transition into a percolation transition of the overloaded grains, and we argue that in the presence of long-range correlations above the transition angle, the granular slope is metastable.Comment: 11 pages, 14 Fig, submitted to PR

The Cohesive Granular Collapse as a Continuum : Parametrization Study

Although intensive research on the flow of dry granular materials has allowed for the proposition of continuum rheology and modelling, the behaviour of flowing cohesive material has attracted less attention so far. To start modelling such cohesive flows, we first focus on the configuration of a granular collapse, which is a simple benchmark test. Specifically, we compare granular-collapse experiments of cohesive grains with numerical simulations, where we test a simple rheology for the material : the so-called µ(I)-rheology, supplmented by a yield stress for cohesion. This document reports the sensitivity of our numerical simulations on the parameters of the rheology, often challenging to measure in experiments

Evolution of nano-pores during annealing of technically pure molybdenum sheet produced from different sintered formats

Molybdenum is a refractory metal with no phase transformation in the solid state and a high melting point. It is therefore an excellent structural material for various high temperature applications. Especially in this field of operation, significant creep resistance is essential. To achieve this, a microstructure with grains in the range of millimeters is desired. However, as demonstrated in the present study, the onset temperature for secondary recrystallization, which would lead to a beneficial grain size, is among other things dependent on the initial dimensions of the sintered part. One possible reason for the different microstructural evolutions is the influence of residual pores in sub-micrometer size. Sheets were thus fabricated via three different production routes employing the same initial Mo powder to exclude chemical variation as an influencing factor. The samples were investigated by in-situ small-angle X-ray scattering at a synchrotron radiation source with two different heating rates. Additionally, selected annealed samples were studied ex-situ with high energy X-rays. The apparent volume fraction of pores is compared to a volatilization model for the vaporization of typical accompanying elements and the induced thermal expansion

Experimental evidence of ageing and slow restoration of the weak-contact configuration in tilted 3D granular packings

Granular packings slowly driven towards their instability threshold are studied using a digital imaging technique as well as a nonlinear acoustic method. The former method allows us to study grain rearrangements on the surface during the tilting and the latter enables to selectively probe the modifications of the weak-contact fraction in the material bulk. Gradual ageing of both the surface activity and the weak-contact reconfigurations is observed as a result of repeated tilt cycles up to a given angle smaller than the angle of avalanche. For an aged configuration reached after several consecutive tilt cycles, abrupt resumption of the on-surface activity and of the weak-contact rearrangements occurs when the packing is subsequently inclined beyond the previous maximal tilting angle. This behavior is compared with literature results from numerical simulations of inclined 2D packings. It is also found that the aged weak-contact configurations exhibit spontaneous restoration towards the initial state if the packing remains at rest for tens of minutes. When the packing is titled forth and back between zero and near-critical angles, instead of ageing, the weak-contact configuration exhibits "internal weak-contact avalanches" in the vicinity of both the near-critical and zero angles. By contrast, the stronger-contact skeleton remains stable

Pre-avalanche instabilities in a granular pile

We investigate numerically the transition between static equilibrium and dynamic surface flow of a 2D cohesionless granular system driven by a continuous gravity loading. This transition is characterized by intermittent local dynamic rearrangements and can be described by an order parameter defined as the density of critical contacts, e.g. contacts where the friction is fully mobilized. Analysis of the spatial correlations of critical contacts shows the occurence of ``fluidized'' clusters which exhibit a power-law divergence in size at the approach of the stability limit. The results are compatible with recent models that describe the granular system during the static/dynamic transition as a multi-phase system.Comment: 9 pages, 6 figures, submitted to Phys. Rev. Let

Memory of the Unjamming Transition during Cyclic Tiltings of a Granular Pile

Discrete numerical simulations are performed to study the evolution of the micro-structure and the response of a granular packing during successive loading-unloading cycles, consisting of quasi-static rotations in the gravity field between opposite inclination angles. We show that internal variables, e.g., stress and fabric of the pile, exhibit hysteresis during these cycles due to the exploration of different metastable configurations. Interestingly, the hysteretic behaviour of the pile strongly depends on the maximal inclination of the cycles, giving evidence of the irreversible modifications of the pile state occurring close to the unjamming transition. More specifically, we show that for cycles with maximal inclination larger than the repose angle, the weak contact network carries the memory of the unjamming transition. These results demonstrate the relevance of a two-phases description -strong and weak contact networks- for a granular system, as soon as it has approached the unjamming transition.Comment: 13 pages, 15 figures, soumis \`{a} Phys. Rev.